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Abstract 

It is predicted that nuclear spin conversion in molecules can be efficiently con- 
trolled by strong laser radiation resonant to rovibrational molecular transition. 
The phenomenon can be used for substantial enrichment of spin isomers, or 
for detection of very weak (10 — 100 Hz) interactions in molecules. 
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I. INTRODUCTION 



Many symmetrical molecules exist in Nature only in the form of nuclear spin isomers 
[0]. These isomers differ by symmetry of nuclear spin wave function and, consequently, by 
symmetry of molecular spatial wave function. Relaxation between different spin states (spin 
conversion) is extremely slow process if it is compared with other gas kinetic rates, e.g., 
vibrational relaxation. It makes nuclear spin isomers unique objects with many potential 
applications. 

It was shown in a recent paper [Q] that external radiation can influence conversion of spin 
isomers both through level populations and through optically induced coherences. Purpose of 
the paper p[ was to investigate main features of the coherent control of isomer enrichment 
and conversion. In order to achieve the goal, the process was considered in a simplest 
arrangement in which microwave radiation excited molecular rotational transition. 

In the present paper we will consider the process in more complicated arrangement which 
promises significantly better control of spin conversion. One outcome of this high efficiency 
is that the coherent control in new arrangement can be used for detection of very weak 
(10 — 100 Hz) interactions in molecules. 

II. EQUATION OF CHANGE 

First, we give qualitative picture of the process. Let us assume that a test molecule 
has two nuclear spin states, ortho and para, and that there is a laser radiation resonant to 
the rovibrational transition, m — n, in the ortho subspace. Fig. |1]. The low state n is not 
mixed with para states, but the upper state m is mixed by the intramolecular perturbation 
V with the para state k. In addition, there is the ortho-para level pair, m' — k', in the 
ground vibrational state mixed by another intramolecular perturbation V'. We assume that 
collisions of the test molecules cannot alter their spin state. This arrangement corresponds 
to the general formulation of quantum relaxation in which one has collisionally isolated 
subspaces of states mixed together by internal perturbation 0. Suppose that the molecule 
is placed initially in the ortho subspace of the ground vibrational state. Due to collisions 
the test molecule will undergo fast rotational relaxation inside the ortho subspace. This 
will proceed until molecule jumps to the state m', which is directly mixed with the para 
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state k', or to the state n which is mixed by combined action of the external field and 
intramolecular perturbation V with the para state k. Admixture of a para state implies 
that the next collision can move the molecule to another para states and thus localize it 
inside the para subspace. It is clear that resonant radiation can significantly modify the spin 
conversion. One knows from the literature on resonant interaction of strong laser radiation 
with matter (see, e.g., @J^) that radiation can change population of states, split states and 
create coherences in the system. All these can affect the mixing of ortho and para states 
and thus the spin conversion process. 

Quantitative description of the problem can be performed with the help of kinetic equa- 
tion for density matrix, p. The molecular Hamiltonian consists of four terms, 

H = Ho + hG + hV + hV'. (1) 

The main part, Hq, has the eigen ortho and para states shown in Fig. |l|. hG describes the 
molecular interactions with the external radiation, 

G = -{Eod/n)cosujLt, (2) 

where Eg and are the amplitude and frequency of the electromagnetic wave; d is the 
operator of the molecular electric dipole moment. We have neglected molecular motion in 
the operator G in order to simplify the theory. 

In the representation of the eigen states of Hq kinetic equation reads, 

dp/dt = {dp/dt),oii -i[G + V + V\p\, (3) 

where {dp/dt)coii is the collision integral. Further, collisions in our system will be described 
by the model standard in the theory of molecular interaction with laser radiation. The 
off-diagonal elements of {dp/dt)coii will be assumed to have only decay terms, 

idpj,j'/dt)cou = -Tpjf, j ^ j'. (4) 

Here, j and j' indicate rovibrational states of the molecule. These states are assumed to 
have no degeneracy. The decoherence rates were taken equal for all off-diagonal elements of 
collision integral. Diagonal terms of the collision integral will be described in the framework 
of the strong collision model. 
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Our goal is to determine time dependence of the total concentration of molecules in one 
spin state. For example, for the total concentration of ortho molecules, Po, one can get from 
Eq. (1^) the following equation of change, 

dpo/dt = 2Rei{pmkVkm + Pm'k'V^'m')- (5) 

In fact, this result is valid for any model of collision integral, as long as collisions do not 
alter the molecular spin state which implies, J2j{dpjj/dt)coii = 0, if j G ortho, or j G para. 

One has to make a few simplifications in order to find the off-diagonal density matrix 
elements, pmk and pm'k'- We assume V and V' being small and consider zero and first order 
terms of the density matrix, 

p = p(°)+p«. (6) 
We start with zero order perturbation theory. p(°^ is determined by the kinetic equation, 

rfpWM=(rfp(o)M),„,-z[G,p(°)]. (7) 
In zero order perturbation theory, para molecules are at equilibrium, 

(^?, j) = in- pi°V(j); (e, j) = 0, (8) 

where n is the total concentration of the test molecules; w{j) is the Boltzmann distribution 
over rotational states. We will assume the same function w{j) for each of four vibrational 
states. We have neglected in Eq. (|^) vibrational resonance exchange between ortho and 
para molecules which would populate the upper vibrational state of para molecules. It was 
taken into account also that ortho-para exchange is on many orders of magnitude slower 
than vibrational and rotational relaxations. 

Equations for stationary level populations of ortho molecules are obtained from Eqs. (^) 
and (|^. Thus one has, 

{lyy + ^R)pf\e,3) = unw{j)p^^\e) + pfp5,^- 

i^RP^^\g.3) = ^RwU)pf{g) + ^vw{3)p^^\e) - p^^^p5,n. 

pT{9,n)-p^^\eM. (9) 



r2 + fi2 

where p^o\^) ^^"^ P^oHd) concentrations of ortho molecules in excited and ground 

vibrational states; z/y and z/^ are the rates of vibrational and rotational relaxations; p is 



the probability of optical excitation of ortlio molecules. Introduction of different relaxation 
rates for different degrees of freedom makes the model of strong collisions more accurate. 
Eqs. correspond to rotational wave approximation. Matrix element of G is given by 

G^n = -Ge-'^*; G = Eod^^/2h, (10) 

where VL = ujl — ^mn is the radiation frequency detuning from the absorption line center, 
oJmn] the line over symbol indicate a time-independent factor. Rabi frequency, G, is assumed 
to be real. 

Solution of Eqs. (P) has no difficulty. Concentration of ortho molecules in excited vibra- 
tional state, p*-°^(e), and in the state m which one needs for further calculations read, 

(0)/ N (0) P (0)/ N (0) P 2TG'^w{n) 

z/;;^ = w{m)uy^ + (1 - w{m)){uv + i^r)'^- (U) 

Here, is the homogeneous linewidth of the absorption spectrum profile, = + 
2FrG^; r = z/"-*^ + u'^ = w{n)h'y^ + (1 — w{n))h'^^. z/^ and z/„ are the effective 

population decay rates of the corresponding states. In a similar way, one can calculate from 
Eq. (|^) the off-diagonal density matrix element, which amplitude is equal to, 

-p^^\m\n)=^Gp^^^w{n)^^. (12) 



In zero order perturbation theory, one neglects perturbations V and V'. It implies that 
re are 
one has. 



there are no coherences between ortho and para states, = 0; p^/^/ = 0. Consequently, 



dpo/dt = 2Rei{p^^lVkm + pSfc'V;',^,), (13) 

instead of Eq. (^. Note, that the spin conversion appears in the second order approximation. 
The first order correction to the density matrix, p^^\ is determined by the equation, 

dp^'^/dt = {dp^'^/dtUi - ^[G,P^'^] - t[V + V',p(°)]. (14) 

For p^)^; one has from this equation, 

&' = ^^[p'^H9,k')-p^^\g,m% (15) 

where u' = oj^'k'- Pmi "^^^ obtained from equations which are deduced from Eq. (0), 



Substitutions, K 



mk 



Ve 



iujt 



rmk 



■tt(I) piut. 
Pmk^ 1 



Pnk 



(16) 



Pnk^^^'^^ '^fc"*)^ transform 



Eqs. (|16]) to algebraic equations from which one finds p^^. Using and p^/fc/ from Eq. ([T5|) 
one has an equation of change ([13|) in the form, 

dpo 2T\V'\^ 



dt 



(r + icu)[r + i(fi + cj)] + G2 



(17) 



III. ENRICHMENT AND CONVERSION 

The denominator of the second term in the right-hand side of Eq. (|17D is convenient to 
present as, (F + iuJi){T + 1002)-, where 



New parameters, Ui and UJ2, can be interpreted as the gaps between the two components, 
I mi > and |m2 >, of the ortho state |m >, spht by the optical field, and the para state \k >. 
The ortho state, \m2 >, crosses the para state, \k >, at Q = —uj{l — G^/u;^), see Fig. 0. 
The splitting of states by resonant laser radiation is well-known phenomenon in nonlinear 
spectroscopy 

Using Eq. (|T^) one can present equation of change in the final form, 

dpo/ dt = n-f'^p - poT^ 7 = + 7^^ - 7n + 7n + Icoh- (19) 

In writing this equation we have neglected in the right-hand side of Eq. ([T7| ) small difference 
between p^^^ and the total concentration of ortho molecules, po- In Eq. (^) the following 
partial conversion rates have been introduced. The field independent rates, 

rl\2 



2F|VP 2V\V' 

7op= p2^^/2 ^(^0; 7;o= (20) 



The rate 7/ree = 7op+7po determines the equilibration rate in the system without an external 
field. The field dependent term. 
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in = ipoVl^V. (21) 

appears because of depletion of the ground vibrational state of ortlio molecules by optical 
excitation, Index n in 7^ comes from "noncoherent", i.e., induced by level populations. 
Another term of similar "noncoherent" origin appears due to the level population, pf^{m\ 
in Eq. (O, 

7,. = 2|>/pAfe .^ + + (22) 
And finally the "coherent" term, 7coh, originated from p^*^-*(n|m), in Eq. (p!7|), 

^^^^ = 2\V\^^Re- -. (23) 

' ' 2r (r + zc^i)(r + 2cu2) ^ ^ 

Solution to Eq. ([T9|) can be presented as, 

Po = Po+ (Po(0) -Po)exp(-7t); Po = ^7op/7- (24) 

Here 7 is the equilibration rate in the system in the presence of external field; is the 
stationary concentration of ortho molecules. Without an external radiation (at the instant 
t = 0), the equilibrium concentration of para molecules is equal to, 

Pp(0) =n- po(0) = ni^J-i^ree, (25) 

if the Boltzmann factors are assumed to be equal, wik') = w{m'). This implies 7op = 7^^ 
(see Eq. (^)), the laser field produces a stationary enrichment of para molecules, 

/3 = ^^-1 = 1-2^. (26) 
Pp(0) 7 ^ ' 

One can see from this equation that external field changes concentration of para isomers if 

7 7^ 7/ree- 

We assume in further analysis the following parameters, u = 100 MHz, uj' = 130 MHz, 
V^,j^, = 5 kHz, r = 2 ■ 10^ s^^/Torr and the Boltzmann factors of the states m', k', m, 
and k all equal 10^^. This set of parameters gives the field free conversion rate, 7/ree = 
10^^ s^^/Torr, which coincides with the conversion rate in ^^CHsF. Nuclear spin conversion 
in these molecules is governed by quantum relaxation (see the review Q). The rotational 
and vibrational relaxation rates will be taken equal, = O.IF and z/y = O.OIF, respectively. 

First, we consider relatively low optical fields, thus small G. In this case one has two 
peaks in enrichment at frequencies Q ^ —00 and f2 = 0, see Fig. 13. The data shown in this 



figure correspond to Vmk = 3 kHz, and F = 2 MHz. The peak at f2 = appears because 
the excitation probability, p, has maximum at this frequency. Amplitude of this peak is 
determined mainly by the rate 7„. As G grows, the amplitude of the peak 2 reaches the 
value 'jn/lfree ~ (Vlj' /2V'ujY which constitutes ~ 15%. The peak at = in isomer 
enrichment was predicted in by considering only the level population effects. 

Peak at ~ —uj appears because the ortho state \m2 > crosses the para state |fc > at 
this frequency of the external field (see Fig. This peak is determined mainly by jcoh- 
When G increases its amplitude grows to much bigger values than the amplitude of the peak 
at = 0. At resonant frequency Q ~ —u the rate ■jcoh is enhanced by large factor (cu/F)^. 
This explains much larger enrichment at i7 ^ —u. Note, that large enrichment occurs only 
if the excitation probability at this frequency, p{—uj), is not very low. 

The data shown in Fig. |^ correspond to strong optical field, G = 50 MHz, and three 
values of Vmk- T was taken equal 2 MHz. One can see, that strong optical field is able 
to convert almost all molecules to the para state if Vmk — V^'k'- Thus relatively weak 
(3 kHz) coupling in upper state is able to produce macroscopic effect, viz., almost complete 
enrichment of spin isomers. It is of fundamental importance, that even for much weaker 
coupling in upper state, enrichment is still significant. For example, if the perturbation in 
upper state, Vmk = 30 Hz, one has the enrichment, (3 ~ 1%. Enrichment at this level can 
easily be measured. It is important that the enrichment peak at 1] ^ —u) is narrow (the 
width ^ F) and thus can be distinguished from much wider structures (the width ~ F^) 
induced by population effects. 

Equilibration rate in the system is given by 7, see Eq. (p4D. It is convenient to characterize 
the conversion rate in relative units, 

Irel = free " 1, (27) 

Conversion rate, like enrichment, has two peaks in its frequency dependence at low G. If 
Rabi frequency, G, is large and the ortho-para couplings in upper and low states have the 
same order of magnitude, conversion can be significantly enhanced (Fig. ^, upper panel). 
Again, this enhancement appears because of the crossing of ortho and para states in upper 
vibrational state by external field. 
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IV. DISCUSSION 



The phenomenon considered in the paper is based on the level splitting produced by 
resonant electromagnetic radiation. Sometimes, this splitting is called in optics the dynamic 
Stark effect. The essence of the effect can be understood as follows. Mixing of ortho and 
para states depends on magnitude of the perturbation V but also on the ortho-para level 
gap. Optical field splits the molecular state and thus change the gaps between the ortho and 
para states. Conversion rate is significantly enhanced when the ortho and para states cross. 
Similar enhancement occurs when ortho and para states are crossed by ordinary Stark effect 
in an external DC electric field [Q. 

In the same way, one can understand high sensitivity of the phenomenon to weak ortho- 
para couplings in excited state. Equilibrium concentrations of ortho and para molecules are 
achieved when one has the ortho-to-para flux in excited state equal to the back flux in the 
ground vibrational state. Back flux is slow because it is determined by the non-degenerate 
ortho-para level pair. On the other hand, the flux in excited state can be significantly 
enhanced by proper choice of radiation parameters which allows to cross the ortho and para 
states in upper vibrational state. 

High efficiency of the proposed enrichment method can be used to detect weak pertur- 
bations in excited vibrational state. First of all, it can be the hyperfine perturbations of 
the same origin and similar magnitude as the perturbations in the ground vibrational state. 
We have seen that hyperfine coupling of the order of ~ 10^ Hz is able to convert almost 
all molecules in one spin state. One can also detect much weaker interactions in molecules. 
In this case one should select the ortho-para level pair in excited vibrational state which is 
not mixed by "ordinary" hyperfine interactions in order to avoid the weak interaction to be 
hidden by stronger, ordinary hyperfine interactions. An interesting case is the crossings of 
states having opposite parity. Mixing of such ortho and para states can be performed only 
by spin-dependent, parity-odd interactions which are not observed in molecules yet. 

V. CONCLUSIONS 

We have performed analysis of the spin isomer enrichment and conversion governed by 
molecular rovibrational excitation. This analysis was done using a few simplifications. We 
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have neglected the Doppler broadening of the absorbing transition, degeneracy of molecular 
states, and resonant vibrational exchange between excited and unexcited molecules. These 
simplifications are not crucial for the existence of the phenomenon. More detailed analysis 
will be done elsewhere. 

We have shown that coherent control of nuclear spin conversion in molecules can be 
efficiently performed by strong radiation resonant to rovibrational molecular transition. A 
possible applications of this phenomenon is the enrichment of molecular spin isomers. An- 
other application is the detection of very weak (10 — 100 Hz) interactions in molecules, which 
can be, e.g., parity-odd interactions. 

ACKNOWLEDGMENTS 

This work was made possible by financial support from the Russian Foundation for Basic 
Research (RFBR), grant No. 98-03-33124a. 



10 



REFERENCES 

[1] L. D. Landau and E. M. Lifshitz, Quantum Mechanics, 3rd ed. (Pergamon Press, Oxford, 
1981). 

[2] P. L. Chapovsky, |http://arXiv.org/abs/physics/00110T2 . 
[3] P. L. Chapovsky, Physica A (Amsterdam) 233, 441 (1996). 

[4] S. G. Rautian, G. I. Smirnov, and A. M. Shalagin, Nonlinear resonances in atom and 
molecular spectra (Nauka, Siberian Branch, Novosibirsk, Russia, 1979), p. 310. 

[5] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, Atom-Photon Interactions (Wi- 
ley, New- York, 1992). 

[6] P. L. Chapovsky and L. J. F. Hermans, Annu. Rev. Phys. Chem. 50, 315 (1999). 

[7] L. V. Il'ichov, L. J. F. Hermans, A. M. Shalagin, and P. L. Chapovsky, Chem. Phys. 
Lett. 297, 439 (1998). 

[8] B. Nagels, N. Galas, D. A. Roozemond, L. J. F. Hermans, and P. L. Chapovsky, Phys. 
Rev. Lett. 77, 4732 (1996). 



11 



m 



n 

m' 



V 



mk 



optical 
excitation 



Vmlc' 



■< ► 



k' 



ortho 



para 



FIG. 1. Level scheme. Horizontal lines indicate the ortho-para mixing in the ground and 
excited vibrational states. There are rotational relaxation inside each vibrational state (rate vr) 
and vibrational relaxation from upper states (rate vy)- 
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FIG. 3. Frequency dependence of the enrichment of para molecules, (3, at G = 1 and 2 MHz. 
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FIG. 4. Enrichment of para molecules, (3, and conversion rate, ^jrei, for three values ofV^k- 
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